A limit on the electron electric dipole moment using paramagnetic ferroelectric Eu 5 Ba 5T1O3 
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We report on the results of a search for the electron electric dipole moment d e using paramagnetic ferroelec- 
tric Eu 5 Bao.5Ti03. The ferroelectric polarization of the crystal creates an effective electric field that makes it 
energetically favorable for the spins of the seven unpaired 4/ electrons of the Eu 2+ to orient along this polariza- 
tion, provided that d e t 0. This interaction gives rise to a magnetization and is therefore equivalent to a linear 
magnetoelectric effect. The ferroelectric polarization present in Eu .5Ba 5 TiO3 increases the size of this inter- 
action compared to other solid-state systems. We obtain d e = (-1.07 ± 3.06 sfcit ± 1.74 sys ) x 10~ 25 ecm, implying 
an upper limit of \d e \ < 6.05 x 10~ 25 ecm (90% confidence). 

PACS numbers: 11.30.Er, 75.85. +t, 32.10.Dk, 14.60.Cd 



The permanent electron electric dipole moment (eEDM) 
has been of experimental interest for nearly half a century 
because it provides a probe of charge-parity (CP) symmetry 
violation in the universe. Through the CPT theorem 0]], the 
existence of a permanent electric dipole moment, which vio- 
lates time-reversal (T) symmetry, would imply violation of CP 
in order that combined operations of CPT are conserved. CP 
symmetry violation is required in the early universe in order 
to explain the currently observed matter-antimatter asymme- 
try I12[]; furthermore, the CP violation in the standard model 
(SM) is not sufficient to explain this asymmetry 01 . Many 
theories that go beyond the SM contain more CP violation 
and therefore predict a larger eEDM that may be detected by 
the next generation of experiments 14J] . 

The traditional method to search for an eEDM involves 
observing precession of an atom or molecule with unpaired 
electron spins in the presence of both magnetic and elec- 
tric fields 0]. This method has been used extensively 0,01 
and has set the best current upper limit on the eEDM of 
\d e \ < 1.08 x 10~ 27 ecm |8-]. Another measurement pro- 
cedure, first suggested by Shapiro [9], involves placing un- 
paired election spins bound to a crystal lattice in an electric 
field. If d e + 0, the electrons will orient along the elec- 
tric field and produce a magnetization in the sample ifioll . 
To date, two experiments produced eEDM limits using this 
approach ll lL 1 1 2fl . The reverse experiment, where the sam- 
ple is magnetized and a correlated polarization is measured, 
has also been performed 1 1311 . These solid-state -based exper- 
iments sacrifice the narrow atomic and molecular transition 
linewidths for a significantly larger signal due to the high den- 
sity of spins present in a solid. 

Perhaps the most important choice for a solid-state eEDM 
experiment is the material. In Refs. lfl4l[l5ll . the advantages 
of Euo.5Bao5Ti03 are detailed over other materials, and a 
short review will be presented here. Euo.sBao.sTiOa has a per- 
ovskite crystal structure and is ferroelectric below 170 K. Our 
samples, which have approximately 65% ceramic density and 
were made in an identical way to those in Ref. lfl5ll . can be 
partially polarized using moderate voltage (< 3 kV or approx- 
imately 20 kV/cm). The magnetic Eu 2+ ions are responsible 
for paramagnetic behavior above approximately 1 .9 K and be- 
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FIG. 1. (Color online) A cut-through schematic of the eEDM exper- 
iment. Note the coordinate system in the bottom of the figure. 



havior consistent with anti-ferromagnetism at lower tempera- 
tures iH . 

The sample magnetization induced by the eEDM is given 

by 



M = 



)( m d e E 



(1) 



where Xm is the magnetic susceptibility, d e is the eEDM of 
the electron, fi a is the magnetic moment of the Eu 2+ ion, and 
E* is the effective electric field. As shown for a similar ma- 
terial, GdsGasOn, the effective electric field is proportional 
to the displacement of the Eu 2+ with respect to the center of 



2 



the oxygen octahedron around it flitl . This displacement has 
been computed to be equal to half that of the displacement 
of the Ti 4+ ions with respect to the 2 ~ IU4I1 and is therefore 
proportional to the polarization of the sample, i.e., E* = kP. 
Using the results in Ref. 11611 and scaling for Euo.sBao .5Ti03, 
we conservatively predict k as (10 MV/cm)/(l //C/cm 2 ). The 
EDM interaction [Eq. [T) can be viewed as a first order, lin- 
ear magnetoelectric (ME) effect in the sample. In this picture, 
the free energy of the sample is modified by a linear term 
a' HP, where a' = Xmd e k/fi a an d H is the applied magnetic 
field. Because the sample is cooled in a zero electric field and 
the experiment is operated at 4.2 K where the sample is para- 
magnetic, both parity and time symmetries are conserved in 
the crystal. A non-zero a' can therefore only arise because of 
the eEDM|0. 

A cut-through schematic of the experimental apparatus 
is shown in Fig. Q] Two disc-shaped samples of diameter 
12.6 mm and height 1.7 mm are held onto a centrally located 
ground plane by two electrodes. Like most of the cryogenic 
components, the ground plane is constructed from G10 fiber- 
glass but is coated with graphite to make the surface conduc- 
tive. An 8-turn superconducting Nb-Ti alloy pickup loop is 
wound inside the ground plane. The pickup loop transfers the 
flux generated by the magnetization of the samples to a super- 
conducting quantum interference device (SQUID) that is used 
as a magnetometer. Because of the geometry of the samples, 
there are demagnetizing fields that lead to suppression of the 
magnetic flux detected by the SQUID 1 1 811 . To electrically po- 
larize the samples, voltage is generated by a custom-built high 
voltage supply and applied via graphite-painted electrodes on 
the flat surfaces of the samples. Additional leads from the 
ground planes attach to a high dynamic range, transimpedance 
amplifier Il9tl . with which currents that flow through the sam- 
ple are measured. The polarization is determined by numeri- 
cally integrating the measured current. Such numerical inte- 
gration is accurate only to an arbitrary constant and thus mea- 
sures the change in polarization but not the absolute polariza- 
tion. 

Two layers of superconducting magnetic shields made of 
1 mm thick, 99.9% pure Pb foil surround the sample region. 
This shielding offers a minimum shielding factor of 10 8 for 
time- varying magnetic fields. However, during cooling of the 
experiment, they trap ambient magnetic fields as they undergo 
the superconducting transition. This trapped field can be can- 
celed using superconducting magnetic field coils wound on 
a cylindrical form of radius 5.2 cm and length 17 cm. A 
solenoid coil applies a field H z parallel to the normal vector of 
the pickup loop (defined to be the z direction), and a cosine-0 
type coil applies a field H x perpendicular to the normal vector 
of the pickup loop but at a random azimuthal angle (defined 
to be the x direction). Lastly, an anti-Helmholtz coil applies a 
magnetic field gradient dH z /dz- 

Fig. |2] shows the experimental measurement procedure. 
Electric field pulses separated by a time r are applied to either 
the top sample, bottom sample, or both to modulate the rema- 
nent polarization. Because Eq. Q] is linear in P, the eEDM- 
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FIG. 2. (Color online) The procedure for measuring the eEDM. Elec- 
tric field pulses (top) of duration t p are applied a time r apart. Each 
subsequent pulse reverses the polarization of the sample. The cur- 
rent flow through the samples (second from the top) is numerically 
integrated to obtain the polarization of the samples (second from the 
bottom). The SQUID signal (bottom) is averaged after each pulse 
between times t s and t, where t s is generally 0.8r. Shown on the 
right are typical orders of magnitude for the various applied fields 
and measured quantities. If d e a 10~ 27 ecm, the size of the SQUID 
signal would be of order 1 n<I>o. 



induced magnetization will be similarly modulated. To mea- 
sure the resulting modulation, the SQUID signal is averaged 
after allowing time for transients to settle. To prevent back- 
ground drifts in the signal from impacting the computation of 
the correlation, the average SQUID signals for four adjacent 
pulses in time are weighted by |, -|, |, -| and summed. 
This procedure determines the difference in the SQUID signal 
between the two polarization states AO and eliminates the ef- 
fect of a linear drift. AO is then divided by the difference in 
the polarization AP to determine the correlation between the 
SQUID signal and the polarization. This correlation AO/AP 
is proportional to the ME coefficient a' and thus d e . 

The predominant noise source is the SQUID magnetome- 
ter's intrinsic noise. Above 1 Hz, the noise spectral density is 
approximately white at 3 pOo/ VMz. Below 1 Hz, the noise 
of the SQUID rises roughly as l/f, where / is the frequency. 
Due to technical constraints, the fastest r corresponds to a re- 
versal frequency of 0.25 Hz, within the l/f noise regime of 
our SQUIDs. Despite operating in the l/f regime of the noise, 
the statistics for a data run are Gaussian. Each data run com- 
prises between 200 and 600 electric field pulses, and a Gaus- 
sian is fit to the distribution of AO/AP. The error of the best 
fit mean is used as the statistical error for that run. The typical 
reduced x 1 f° r suc h a fit is near unity. Because the samples 
are reversed at a frequency within the l/f noise regime of the 



3 



e 



|20 

a 

15 



10 



1 

J' yr^S- 


i 
i 
i 


i 




f 













20 
10 


-10 

10 



-10 
-20 



2 



2 4 6 

Time (s) 

FIG. 3. (Color online) Example of the difference of the heating decay 
transient. The blue dashed lines show the applied electric field pulses 
and the red solid lines show the resulting SQUID signal. The large 
features in the SQUID signal seen during the electric field pulses 
are caused by the current that flows during the polarization reversal. 
After the reversal, the heated sample returns to equilibrium with the 
LHe bath, which can be seen as the decay after the pulse. These data 
were taken in the presence of an H x field (top panel) and an H z field 
(bottom panel), each approximately 1 mG. 



SQUIDs, the statistical errors of A<t>/AP tend to be an order 
of magnitude larger than those projected in Ref. 11511 . 

Several systematic effects in the experiment can generate a 
non-zero A<t>/AP and can therefore mask or mimic the linear 
ME effect due to the eEDM. For example, if the samples are 
in a non-zero magnetic field, a change in the temperature of 
the sample(s) will lead to a change in permeability that will 
subsequently change the flux through the SQUID. Because of 
the dissipation inherent to ferroelectrics, polarization reversals 
heat the sample(s). As the samples return to equilibrium with 
the liquid helium bath, a transient can be seen in the SQUID 
signal, as shown in Fig. [3] Provided this heating is equal when 
the sample polarization is switched from +z to —z (a negative 
pulse) and —z to +z (a positive pulse), the heating transients 
are identical for positive and negative remanent polarizations, 
and there is no systematic effect. A measure of the amount of 
heat released by a given pulse can be derived from the integral 
of P ■ dE, where P is the polarization and E is the applied 
electric field lT20ll . and is of the order of 1 mJ per pulse. 

To quantify the size of the resulting A<J>/AP, magnetic fields 
were applied and the electric field pulses were deliberately 
unbalanced to produce different heating for positive and neg- 
ative pulses. The resulting correlation was measured in this 
manner for each reversal frequency and for each sample. The 
correlations were then fit to AQ/AP = aAQp, where AQ is 
the difference in heat released between a positive and negative 
pulse, p is a proxy for the magnetic field, and a is a tunable 
constant. As shown in Fig. [3] the transient is significantly dif- 
ferent for H x vs. H z fields; for this reason, the fits for the 
correlation use p = (d<$>/dt) as a proxy for the strength of H z 



and p = (d 2 <$>/dt 2 ) as a proxy for the strength of H x . The re- 
sulting fits to experimental data confirm the validity of these 
proxies. The best fit values for a are used to predict the size 
of the correlation when the magnetic field is close to zero and 
the electric field pulses are symmetric. In this configuration, it 
is not known a priori what type of field envelops the samples; 
therefore, the most likely correlation for both an H z field and 
an H x field is computed. The resulting predictions are used as 
a 1-cr systematic error without applying any correction. 

In addition to this heating effect, the higher-order ME ef- 
fect that is present in titanates can also generate a non-zero 
A<t>/AP. Given the symmetries present in our sample, the 
magnetization induced by the higher-order ME effect will be 
given by M — 5x m P 2 H, where P is the absolute polarization. 
Using the same experimental apparatus, the constant 6 was 
measured for this material; details will be presented in a later 
paper in preparation. Because the magnetoelectric-induced 
magnetization depends on P 2 , a non-zero correlation will re- 
sult only if the two different absolute polarization states in the 
modulation have different magnitudes. Thus, the error in de- 
termining the absolute zero of polarization will determine the 
maximum possible difference in P 2 when the polarization is 
reversed. The error in the absolute zero of P is taken to be 
0.1 yt/C/cm 2 at 95% confidence, which is motivated by the fi- 
delity with which samples can be depolarized using electric 
fields. Depolarization effectively resets the constant of inte- 
gration in the determination of the polarization and thus the 
fidelity limits our knowledge of the absolute zero of the po- 
larization. Using this error estimate for the absolute measure- 
ment of P, a A<t>/AP is computed and used as a systematic 
error. 

Because of the inherent dissipation present in ferroelectrics, 
the sample takes some time to reach the final polarization state 
after the electric field is applied. This phenomenon is known 
as dielectric relaxation 112 If! . As the sample relaxes to its fi- 
nal state, current continues to flow through the sample. This 
current scales as f _1 , where t is the time since the polariza- 
tion reversal. To suppress this dielectric relaxation, an addi- 
tional time-varying voltage (maximum 40 V) is applied using 
a proportional-integrator-differentiator (PID) circuit to force 
the net current to zero. To estimate a AcP/AP that may result, 
the SQUID response during the electric field pulse is used to 
calculate the sensitivity of the SQUID to the current through 
each sample. The effect on the SQUID signal due to any cur- 
rent that is not suppressed by the PID is then computed and 
used to estimate the correlation. The correlation due to di- 
electric relaxation is then taken to be a 1-cr systematic error. 

The total integrated time for the data used in the final analy- 
sis is approximately 1 hour and 40 minutes. All data where the 
same sample(s) are driven at the same reversal frequency and 
with the same amplitude electric field pulses were averaged 
together, weighted by their statistical errors. The correlation 
is then converted into a linear ME coefficient and an equiv- 
alent d e . To enable comparision with linear ME coefficients 
that are expressed in units of s irT 1 , we define a = Xe^a' ', 
where Xe = PI £qE ~ 700 is an effective electrical susceptibil- 
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FIG. 4. (Color online) Final best fit eEDM values with statistical 
error bars. Data were recorded at reversal frequencies of 0.25 Hz 
(violet, upward triangles), 0.167 Hz (red squares), 0.125 Hz (green, 
downward triangles), and 0.1 Hz (blue circles). The annotations 
show which sample(s) were driven. The solid cyan line shows the 
best fit mean, and the dashed cyan lines show its l-a statistical error. 
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Heating (H x ) 


0.71 


1.82 


-0.07 


Heating (H-) 


-0.05 


-0.72 


-0.02 


Dielectric relaxation 


-0.97 


-0.11 


0.70 


Higher-order ME effect 


1.40 


0.26 


0.47 



TABLE I. Breakdown of the systematic errors in the experiment 
by source and which sample(s) were driven. Units of the table are 
10~ 25 ecm. 



ity for the ferroelectric. The results are shown in Fig. @] and 
are consistent with zero. 

A breakdown of the systematics is shown in Tab. U The 
higher-order ME effect produces a significant systematic ef- 
fect because of the conservative estimate of our knowledge of 
the absolute polarization. The systematic due to the heating 
shows complicated behavior, and is significantly less when 
both samples are driven. The reason for this reduction is 
twofold. When driving both samples, there is a significant 
rejection of the effect of a transverse field because H x couples 
to the two samples in the opposite way. Second, the asymme- 
try in the heating when both samples were used was measured 
to be nearly equal and opposite, leading to rejection of H z . 

The final best fit results are d e = (- 1 .07 ±3.06 sta t± 1 -74 svs ) x 
1(T 25 ecm and a = (-0.57 ± 1.64 stat ±0.93 sys )X 1(T 21 s/m Q. 
This result implies an eEDM limit of \d e \ < 6.05 x 10~ 25 ecm 
(90% confidence). Compared to previous solid state eEDM 
measurements, this limit is approximately a factor of ten im- 
provement over Ref. 1 1311 and a factor of three better than 
Ref. d|. 

In conclusion, we have built and operated an experiment 
that has established an upper limit on the eEDM better than 
any previously published solid-state experiment. The use of 
Euo sBao 5TiC»3 takes advantage of the ferroelectric polariza- 
tion of the crystal to enhance the effective electric field inter- 
acting with the EDM, which is similar to the mechanism that 



enhances the sensitivity of diatomic polar molecule searches. 
The ultimate EDM limit can be improved in future versions of 
the experiment by identifying and suppressing the sources of 
excess noise in the SQUID magnetometers below 1 Hz. Fur- 
ther suppression of systematics, such as heating and dielectric 
relaxation, could be obtained by improving magnetic shield- 
ing and optimizing the current feedback system. Alternatively, 
these systematics may be suppressed by using a low-loss para- 
electric doped with Eu 2+ ions, such as SrTiC<3 M23L 12411 . 
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